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Abstract 

Ultra High Energy Cosmic Rays with energies above ~ 10 18 eV provide an unique window to study 
hadronic interactions at energies well above those achieved in the largest man-made accelerators. 
We argue that at those energies string percolation may occur and play an important role on the de- 
scription of the induced Extensive Air Showers by enhancing strangeness and baryon production. 
^ This leads to a significant increase of the muon content of the cascade in agreement with recent 
Q-i data collected at UHECR experiments. In this work, the effects of string percolation in hadronic 
interactions are implemented in an EAS code and their impact on several shower observables is 
evaluated and discussed. 

Keywords: Extensive Air Shower, Muon content, Strangeness Enhancement, String Percolation 
Model 



1. Introduction 

Ultra High Energy Cosmic Rays (UHECRs - with energies above ~ 10 18 eV), are the most 
energetic particles known in the Universe yj]. These rare subatomic particles reach the Earth 
and interact with the nuclei in the atmosphere initiating huge cascades of lower energy particles 
usually referred in the literature as Extensive Air Showers (EAS). The origin and composition of 
the UHECRs are still unknown and are an intense research subject Their interactions in the 
atmosphere are also very interesting, as the centre-of-mass energy of the collisions reaches some 
hundreds of TeV, well above the energies reached at the largest man-made accelerators to study 
particle physics. This unique physics window can be indirectly explored in large arrays devoted 
to UHECR detection, of which the Pierre Auger Observatory at Malargiie, Argentina yfl is the 
paradigmatic example. 

Recent results on the muon content of EAS at ground level obtained at the Pierre Auger Ob- 
servatory, have demonstrated that the simulations of EAS show a significant deficit in the number 
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of muons predicted when compared to data [|4J, |5fl . The relative deficit is ~50% when data are 
compared to simulations of proton-induced showers with the QGSJet-II model. This discrepancy 
cannot be explained by the unknown UHECR composition alone. On the other hand, with data 
collected at the KAS CADE- Grande experiment |@] it has been shown that the muon content at 
ground level up to energies of ~ 10 17 7 eV is, as expected, between the number of muons predicted 
in proton and iron initiated showers, regardless of the high energy interaction model being used. 
This scenario would necessarily imply a transition in the behaviour of the hadronic interactions 
in the energy region between the highest energy events in the KAS CADE- Grande array (~ 10 17 7 
eV) and the lowest energy of the events used for muon analysis at the Pierre Auger Observatory 
(~ 10 18 - 4 eV). 

The increase of the number of muons in simulations of EAS can be achieved by effectively 
reducing the amount of shower energy that is transferred to the electromagnetic component, which 
is mainly attained by the decay of neutral pions 0]. 

The string percolation model (SPM) has been quite successful in describing global features of 
data in AA collisions at RHIC (with ene rgy u p to 200 GeV) and in pp collisions at LHC (with 
energy up to 7 TeV for the moment) 11,19, 10]. In the model QCD strings are formed, stretching 



between beam and target. The strings may interact by fusion of clusters of overlapping strings in 
the transverse plane of the collision, in a situation similar to two-dimensional percolation theory 



II11LI12L I13I1. The predictions of the SPM are: i) a reduction of the height of the rapidity plateau 
around y = 0; ii) an increase in the transverse momentum p T of the particles produced in the 
collisions due to string interactions; iii) heavy quark production as a consequence of clustering 
and percolation. In this paper we argue that string percolation effects may be found at very high 
energies ( s/s ~ 40 TeV) in UHECR data. Percolation effectively results in an enhancement of kaon 
and baryon production which results in the reduction of the amount of energy being transferred 
to the electromagnetic cascade in the high energy interactions. It is important to note that the 
percolation model does not forbid some other phenomena such as a rapid increase of the primary- 



air cross section, as for instance described in lll4h . This would actually drive the behaviour of the 



momenta of the distribution of the depth of maximum of the electromagnetic component (Xjj^) of 
the shower. 

The paper is organized in the following way: in Section |2] the string percolation model is 
presented and its impact on the hadronic interaction physics is discussed; afterwards (in Section[3]) 
we explain how the percolation effects were implemented in an EAS code; in Section @] the effect 
of percolation on shower observables, in particular the depth of shower maximum and the number 
of muons at ground, is evaluated and discussed. Some final remarks and conclusions are given in 
Section [5] 

2. String Percolation Model 



In the string model (see for instance lll5L ll6n). colour strings stretched between the beam and 



the target are formed in high energy collisions. In the impact parameter (transverse) plane of 
a collision the strings appear as disks with a characteristic transverse radius r . The disks may 
overlap, cluster and fuse (percolate). The relevant parameter in the String Percolation Model 



(SPM) is the so-called transverse string density rj defined as IU7I1 . 
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where R is the effective radius of the interaction area in the transverse plane, and N s is the average 
number of strings formed in the collision, which typically grows with the centrality and energy of 
the collision tl8fl. Depending on the value of 77 the strings may occupy a significant fraction of 
the transverse area of the collision, overlap, cluster and percolate. If rj < 1 the average number 
of strings per cluster (AO is approximately one (the strings do not overlap), while if 77 » 1 then 
(A 7 ) ~ N s (and the majority of the strings overlap). It is then important to relate the average 
number of strings per cluster (AO to the average area (A) occupied by a cluster. Following [19] (A) 
is written as, 

1 



(N) = (A) 
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with (A) given in units of r by, 
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where f{rf) is a percolation function lll8L ll9ll. 

/O7) = (1 +e-^' )/a )- 1 . 



(4) 



f(rj) characterizes the transition from individual to percolated strings, controlled by the parameters 
77c, the critical transverse string density, and a, the parameter controlling the slope of the curve at 
the transition point. Comparison with Monte Carlo simulations Ill9ll gives 7] c = 1.15 and a = 0.07 
in the case of a uniform proton density profile as a function of the impact paramater. For a Gaussian 
profile the parameters would be rj c = 1.59 and a = 0.14. 

One should note that, contrary to A^, the average number of strings per cluster, and the quantity 
(A) are sensitive to the percolation transition. If 77 — > then (A) — » 1 (in units of r ) and when 
rj -> 00 then (A) -> (R/r Q ) 2 . 

In order to relate the average number of formed strings N s to the collision energy -^s, we use 
the parameterisation in Iil8ll and write, 



N s = b + (2-b)\- 

^0 



(5) 



where the free parameters b = 1.37, -\[sq = 10 GeV and A = 0.23 have been adjusted to experi- 



mental data as explained in II18L |2( 

According to the string percolation model, when strings in a collision fuse the length of the 
rapidity plateau Ay increases in proportion to the number of clustering strings. Following the 
treatment in Ill8ll . 

Av~21og((A0) (6) 
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Figure 1 : Length of the rapidity plateau in a high energy collision AY (in rapidity units) as a function of -\fs. The red 
line is the result for a proton with a uniform density profile distribution, while the blue line corresponds to a Gaussian 
profile. The beam rapidity is also shown in green for reference. All the results are displayed only for the positive 
(forward) rapidity region of the collision. 
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In Fig. [H the increase in the length of the rapidity plateau as a function of V* given in Eq. © is 
plotted for the forward region of the collision only. The beam rapidity is also plotted for compari- 
son. 

Since the strings are stretched in rapidity, energy conservation implies a natural reduction of the 
hadronic interaction multiplicity. The most forward and backward regions of the collision remain 
untouched, i.e., the rapidity of the strings there is not increased. The reduction of multiplicity is 
realized through a colour summation factor F(rj) decreasing with 77 and given by lll6L l2lll. 



F(tj) = 



1 -g-" 



(7) 



In the spirit of the Schwinger mechanism for particle production Il22ll the probability of pro- 
duction of a qq pair (meson) in a string is, 



Pqq <* £ 



(8) 



where m q is the corresponding quark mass and k 0.2 GeV 2 is a phenomenological parameter 
related to the effective color field, the so-called string tension 112311 . If we compare the probability 
of producing a heavy quark pair Pqq (such as a strange quark pair for instance), to the production 
of a light quark pair P qq (such as a uu pair for instance), we obtain, in the absence of percolation: 



Psuin = 0) = ^(77 = 0) ~ exp 

* 11 u 



(9) 



When string percolation is effective the string tension increases due to the presence of higher 
color charges associated to the clusters. It is natural to associate the increase of the string tension 
to the increase of the number of strings per cluster corrected by the color reduction factor F{rf) in 
Eq. ©. With percolation we have an increased string tension k given by, 



k = k (N)F(i 1 ) 



(10) 



As a consequence the probability of producing ass pair relative to the probability of producing a 
uu pair is given by, 

--M-ml)} (11) 



— (17) ~ exp 

uu 



which can be cast as, 



—(77) ~ exp 

L UU 
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(12) 



where P su (0) is the probability of producing a quark s relative to a quark u without percolation. In 
the hadronic interaction models this quantity is of the order of 0.1, in particular in the QGSJet-II 
model P su (0) = 0.07. 

Combining Eqs. CEHS]), © and (fT2l the evolution with energy of the probability of producing 
strangeness in QGSJet-II as a function of -y^ can t> e obtained. The ratio P S slPuu is shown in 
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Fig. [21 The string percolation model predicts a transition in the center-of-mass energy region of 
V* ~ 20 - 30 TeV, above LHC energies, and typically below the energies of the events collected 
at the Pierre Auger Observatory (the energy of Auger events used for muon analysis is above 
E = 10 18 4 eV, which in the center of mass of proton-proton collisions corresponds to V* - 68 
TeV). 



E[eV] 




10 2 10 3 10 4 10 5 10 6 

is [GeV] 

Figure 2: Probability of producing a ss relatively to uu pair as a function of ^[s in the QGSJet-II model accounting for 
string percolation effects. The red line is for a proton with an uniform density profile while the blue line is assuming 
a Gaussian profile. The upper scale is the equivalent energy of a proton in the laboratory frame of the Earth. 

Similarly to the production of qq pairs, in the string percolation model the probability of pro- 
ducing a baryon (associated to the diquark production, qqqq) gets increased in essentially the same 
way, due to the increase of the string tension of the fused string. As a consequence an increase in 
the number of baryons is expected. 

The parameters of the model are presented in Table [TJ The two first parameters, A and r/R, 
are given by the experimental accelerator data and the best fit values are displayed in the column 
designated Preferred Value. The parameters allowed region, due to the experimental uncertainties 
are also shown. The remaining two parameters, r\ c and a, are obtained through simulations and the 
presented values reflect the uncertainty on the proton density profile as a function of the impact 
parameter. Two extreme case are taken: Gaussian and uniform distribution. 
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Parameter 


Min Value 


Max Value 


Preferred Value 


A 


0.21 


0.28 


0.23 


r/R 


0.2 


0.3 


0.2 


1c 


1.15 


1.59 




a 


0.07 


0.14 





Table 1: Main parameters of the string percolation model (see text for definitions). 



3. Percolation of strings in Extensive Air Showers 

To study the impact of percolation of strings in high energy collisions on the properties of 
UHECR-induced showers, we used the EAS simulation code CONEX 0, Q. CONEX is a 
hybrid Monte Carlo for the study of EAS. In CONEX the most energetic interactions, those more 
relevant for the shower development, are simulated explicitely with the Monte Carlo, while the 



remaining interactions are treated by solving cascade equations H24U 25H. This makes it a very fast 
code allowing the simulation of large statistical samples of EAS. However, with CONEX only 
information on the shower longitudinal development along shower axis (X max , number of particles 
at a given depth) can be obtained, and there is no information on the lateral spread of the shower 
particles in the direction perpendicular to the shower axis. Despite this, CONEX is the ideal 
tool for studying the effect of percolation on the average behaviour of some of the main shower 
observables, and their dependence on shower energy - a task that naturally increases the amount 
of EAS simulations needed and calls for a fast and flexible simulation of EAS. 

The hadronic interaction model used as a baseline for the study of the influence of percolation 



of strings on EAS development is the QGS Jet- 11.03 model H26L I27H. The percolation mechanism 
was implemented in QGSJet-II.03 following two steps: the probability of heavy quark production 
was enhanced according to Eq. (fT2l) ; and the rapidity of the strings was increased following Eq. ©. 
The latter step leads to a reduction of multiplicity in the collision as will be explained below. All 
these changes were done in a way that does not compromise the consistency of the hadronic 
interaction model. This means that the modifications to QGSJet-II.03 were done either on input 
parameters of the code (the increase of the probability of heavy quark production), or by applying 
ad-hoc changes to the particles produced by the model after a hadronic interaction (as done to 
emulate the behaviour of the increase of string rapidity and the reduction of multiplicity). All 
the hadronic interaction features described by the model remain untouched. Below we give more 
details on how these changes were implemented. 

The probability of producing strange quark pairs with respect to light quarks pair is one of the 
input parameters in the QGSJet-II.03 model. In the original code, this parameter does not depend 
on the energy of the collision and is fixed to P su (0) = 0.07. To implement the effects of percolation 
of strings, the value of P su was modified using Eq. (fl"2l) . As explained before, P m depends on the 
interaction energy and this dependence is plotted in Fig. [21 As argued in Section[2l the increase of 
the string tension should also affect the probability of diquark production (and hence the produc- 
tion of baryons), and in fact we have assumed the same evolution with Vs for the enhancement of 
the probability of diquark pair creation. The impact of such modifications on particle production 
in high energy hadronic interactions (at a reference energy of ~ 100 GeV) can be seen in the 
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left panel of Fig. [3] As expected there is a notorious increase in the number of strange mesons 
being produced, in particular kaons, as well as in the production of baryons with respect to the 
QGSJet-II.03 model without string percolation. Electromagnetic particles and particles responsi- 
ble for the production of electromagnetic particles (mainly neutral pions) are suppressed. This is 
relatively important for the electromagnetic component of the UHECR-induced showers since the 
decay of n°s is the main channel that feeds this shower component. 




Type 

Figure 3: Left panel: Probability of particle production as a function of the type of particle as predicted with the 
QGSJet-II.03 model without percolation ("Normal QGSII") and after enhancing the probability of strange quark 
production as expected in the percolation model and explained in the text. Right panel: Energy spectrum of the 
produced particles as predicted with the QGSJet-II.03 hadronic model without percolation and after increasing the 
rapidity of the strings as expected in the percolation model and explained in the text. In both cases proton-air collisions 
at E — 10 19 eV were simulated. 

The output of each hadronic interaction is obtained with the QGSJet-II.03 model without per- 
colation. The resulting particles are then processed in order to emulate the effects of percolation on 
the increase of string rapidity. For this purpose we have worked in the collision centre-of-mass ref- 
erence frame and particles are then boosted to the laboratory frame and injected into the CONEX 
Monte Carlo code. Only the properties of the particles with rapidity |y| < yb eam - Ay, (where ybeam 
is the beam rapidity and Ay is the rapidity increment in Eq. ©), were modified to implement the 
effects of percolation. This last condition is necessary to fulfill energy conservation, and at the 
same time it leaves the leading particle (and diffractive particles) untouched, where percolation 
effects are not expected to play a role. Focusing the discussion on positive rapidities the average 
value of the rapidity is increased according to Eq. ©. The rapidity distribution is divided in N 
bins and each bin is stretched in rapidity by an amount Ay/N. Afterwards, the particles in each old 
bin in rapidity are randomly distributed in the new bin so that their rapidity increases on average. 
In each bin, local energy conservation is applied: particles in the old bin in y are redistributed into 
the new stretched bin until the total energy of the particles in the old bin is completely spent. This 
leads to a natural reduction of the number of particles in each bin, and hence to a reduction of the 
total multiplicity of the collision. The final result is a smooth energy distribution that, in the center- 
of-mass framework, has a higher average energy but less particles due to energy conservation. The 
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impact of this implementation on hadronic interactions can be seen in the energy spectrum of the 
particles in the laboratory frame shown in the right panel of Fig. [3] Several structures emerge: the 
peak at lower energies is due to the collision spectators (note that these are proton-air collisions), 
while the peak at the highest energies is attributed to diffractive events. Percolation is not expected 
to change these two features, and they are equal before and after the effects of percolation were 
implemented. As the laboratory framework is highly boosted in one direction relative to the centre 
of mass frame, the particles are either migrated to higher or lower energies, depending on their 
direction with respect to the boost. One can also easily see the reduction of multiplicity due to 
percolation. At E = 10 19 eV this reduction is ~ 65%, i.e. the multiplicity of a hadronic interaction 
with percolation would be in average ~ 35% of that in a standard interaction. 



4. Percolation effects on shower observables 

In this section we present and discuss the impact on shower observables of the presence of 
string percolation in hadronic interactions at ultra high energies. 

As pointed out in Section [31 one of the main consequences of string percolation in hadronic 
interactions is a reduction of both electromagnetic particles and of particles such as neutral pions 
that when decaying feed the electromagnetic component of the shower (see left panel of Fig. 0. 
The energy in these particles will be transferred to the hadronic component of the cascade, and a 



small fraction of it will contribute to the so called invisible energy ll28Ll29h . defined as the energy in 



the shower carried by muons and neutrinos. This effect has consequences on the energy calibration 



performed in modern UHECRs experiments, such as the Pierre Auger Observatory [30] and the 



Telescope Array experiment 113 111 . In both experiments, showers observed simultaneously with 
the fluorescence and the surface detector (typically ~ 10 - 15% of the events) serve to calibrate 
the energy of the showers detected with the surface array alone. In fluorescence detectors, the 
fluorescence light emitted by the excitation of the nitrogen molecules due to the development 
of the shower through the atmosphere is collected, providing in this way a gwas7-calorimetric 
measurement of the primary energy. The excitation of the nitrogen molecules is mostly induced by 
low energy electrons. String percolation should affect the energy assignment by the fluorescence 
detector since, as stated above, the energy transferred to the electromagnetic particles in the shower 
decreases when percolation of strings is effective. This effect is shown in Fig. |4] where the ratio of 
the average electromagnetic energy with percolation, ((E) {p) ), and without percolation as obtained 
in Monte Carlo simulations of EAS performed with CONEX, is plotted as a function of the primary 
cosmic ray energy. The ratio of the electromagnetic energy in showers with percolation to standard 
proton showers, simulated with QGSJet-II, is seen to decrease as the primary energy increases. At 
E = 10 19 eV the energy estimated through the longitudinal development of the light collected by 
the fluorescence detector would be ~ 5% smaller when percolation is effective, inducing in this 
way a systematic underestimate of the energy scale used to calibrate the events observed with the 
surface detector alone. At lower energies, E = 10 17 eV, the percolation transition has not occurred 
yet (see Fig. [2]) and the ratio in Fig. [4] is one by definition. In Fig. |4] the results for two extreme 
cases of the proton density profile: uniform distribution and Gaussian distribution are shown. In 
both scenarios the trend is the same. 

The impact of percolation on the number of muons at ground level is shown in Fig. [5] where the 
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Figure 4: Ratio of the average energy transferred to the electromagnetic component of showers with percolation 
((E)(p)) and without percolation effects (normal QGSJet-II) as function of the energy of the primary proton inducing 
the shower. The red solid line is for a proton with a uniform density profile distribution, while the blue dashed line 
assumes a Gaussian density profile. 
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ratio of the number of muons in showers with percolation effects and in standard proton showers is 
plotted as a function of the primary particle energy. As expected, the enhancement on the number 
of strange mesons and baryons, produced in the highest energy collisions in the shower, leads to 
an increase of the average number of muons at ground level. This effect accounts for ~ 30% of 
the increase, while the hardening of the particle energy spectrum due to percolation (right panel of 
Fig.0, which leads to a more penetrating shower, accounts for an additional increase of ~ 10%. 
The total increase on the number of muons at ground level at E = 10 19 eV is ~ 34-43% depending 
on the percolation model assumptions. The effect increases with the primary energy, both for a 
uniform density profile of the proton (left panel of Fig. [5]) and for a Gaussian proton density 
profile (right panel). Therefore, the string percolation mechanism, provides a possible explanation 
for the deficit of muons in standard simulations observed at ultra-high energies E > 10 18 eV when 
compared to experimental data [0]. Moreover, given the energy at which the onset of percolation 
phenomena is expected, it becomes clear why standard hadronic interactions models are able to 
bracket the lower energy data observed for instance in the KAS CADE- Grande experiment, where 
a deficit is not seen because string percolation is not effective. 
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Figure 5: Ratio of the average number of muons in showers with percolation ({N^ ^) and without percolation effects 
(simulated with the standard QGSJet-II) as function of the energy of the primary proton inducing the shower. The 
left panel corresponds to a uniform density profile distribution, while in the right panel a Gaussian density profile is 
assumed. The impact of the energy scale on the ratio is shown with the green dotted line. 



The reduction of the energy transferred to the electromagnetic component would effectively 
lead to a lowering of the energy scale, as explained above. This means that a number of muons 
measured at ground level would be attributed to a slightly lower cosmic ray energy than in reality 
if percolation effects are not accounted for when estimating the energy of the event. This effect 
is shown in Fig. [5] (dashed lines) where a shift of the energy scale following the trend in Fig. [4] is 
applied. As a consequence at E = 10 19 eV the number of muons with respect to standard proton 
QGSJet-II showers is ~ 40 - 49% higher. 

It is important to notice that the increase in Fig. [5] affects the total number of muons arriving 
at ground level, integrated over all distances to the shower axis. Thus, the features observed in 
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Fig. |5] are not expected to be directly comparable to those observed by the Pierre Auger Observa- 
tory where the number of muons at a distance of 1000 m from the core was measured. A direct 
comparison would require a precise knowledge of the shape of the lateral distribution which is 
not possible to obtain in simulations with the CONEX Monte Carlo, and would require a full 3D 
simulation (such as for instance CORSIKA) which is currently out of the scope of this paper. Fi- 
nally, it is also important to remark that we do not intend to argue about the absolute number of 
muons at ground, but only to discuss general trends with energy, as the number of muons still has 
some dependence on several low energy hadronic interaction phenomena which are not yet fully 
understood [321. 



x10° 




Figure 6: Upper panel: evolution of the number of muons with the shower zenith angle 9. Blue open dots: showers 
with percolation. Red filled dots: showers without string percolation. In both cases proton induced showers at 
E = 10 19 eV were simulated. Bottom panel: ratio of the average number of muons at ground level with and without 
percolation. These results were obtained assuming a uniform density profile for the proton. 



The dependence of the total number of muons at ground level with energies above 1 GeV, with 
the shower zenith angle 6, was also investigated and the results are displayed in Fig. [6j In the 
upper panel of Fig. [6] the number of muons as a function of 8 is displayed for standard showers 
(without string percolation) and for showers with string percolation effects included. The number 
of muons is larger in showers with string percolation effects as explained before and seen in Fig.|5l 
but in both cases (percolation and standard) the total number of muons at ground level follows the 
same trend and decreases with 6. In fact, the ratio between the number of muons in showers with 
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percolation and standard showers is constant with 6 as shown in the lower panel of Fig. [6l These 
plots were done for the case of an uniform density profile of the proton but the same conclusion 
can be withdrawn for a Gaussian profile. This means that percolation phenomena alters the total 
number of muons at ground, but its evolution with the shower zenith angle remains the same. 




log(E/ev) log(E/ev) 

Figure 7: Difference between the average depth of maximum of the electromagnetic and muonic longitudinal profiles 
X^™ x and X^a* respectively, as obtained in proton-induced shower simulations with ((x^^)) and without ((x%^\) 
percolation effects implemented. The effect of the percolation energy scale (see text for explanation) is shown with 
dashed lines. Left panel: an uniform density profile for the proton is assumed. Right panel: a Gaussian profile is 
assumed. 



The effect of percolation on the depth of maximum of the electromagnetic and Muon Produc- 
tion Depth profile |33j,|34J], X^™ and X^ ax respectively, was also studied. The results are displayed 
in Fig. |7] as the difference between the average depth of maximum in showers with percolation 
and the corresponding average depth in standard showers, assuming an uniform density profile 
for the proton (left panel) and for a Gaussian profile (right panel). In both cases the difference 
in the depth of maximum can be seen to increase with primary energy when percolation is effec- 
tive for both X^ x and X^ ax . This is easily understood as the string percolation model predicts 
collisions with lower multiplicity and faster particles on average, which transport more energy in 
the forward region of the collison, pushing both X^ x and deeper into the atmosphere than 
in the case of standard showers without percolation effects. At E = 10 19 eV the value of X^ x , 
typically observed in UHECR experiments with the fluorescence detector, is ~ 22 - 31 g cm -2 
deeper in proton-induced showers with percolation effects than in standard proton-induced show- 
ers simulated with QGSJet-II. Although the effect is not very large, it may alter the interpretation 
of UHECR data in terms of composition of the primary particles B35Q- The depth of the muonic 
longitudinal profile at production X^ nax is also increased by ~ 33 - 47 g cm -2 with respect to that 
obtained in standard showers. Interestingly, the difference between Xf nax in showers with and with- 
out percolation increases continuously with the primary energy while the corresponding difference 
for X^ x saturates at the highest energies. This indicates that the hardening of the energy spectra 
and the increase of strange particles and baryons, have a different impact on the electromagnetic 
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and the hadronic shower components. Moreover, notice that the probability of producing a given 
type of particle saturates quickly after the percolation transition (see Fig. El), while the increase of 
the string rapidity due to string fusion, continues to grow steadily (see Fig.Q]). 

The impact of a lower energy scale as predicted in the percolation model is also shown in 
Fig. [71 although the effect is very small. X^ x increases by only ~ 2 g cm" 2 . The small effect of the 
energy scale on the X max , contrary to what happened with the number of muons, is a consequence 
of the logarithmic X max evolution with the primary energy. 

The results on the electromagnetic X^" x and the muonic X^ nax provide a very distinct exper- 
imental signature which could in principle be searched in hybrid experiments such as Pierre 
Auger Observatory or the Telescope Array Experiment. These experiments combine two inde- 
pendent detection techiques to investigate EAS: Fluorescence Detectors and Surface Detectors. 
The Fluorescence Detectors Oo. l35l. Bill, observe the longitudinal shower development, in par- 
ticular X* n ™ x , by detecting the light emmited by the nitrogen molecules as the shower secondary 
particles (mainly electrons) travel through the atmosphere. The surface detector sample the shower 
secondary charged particles that arrive at ground. Amongst these particles are the muons, which 
can be used, along with their timing information, to reconstruct the Muon Production Depth [37], 
giving access to X„ ax . 

Note that in the percolation scenario the evolution of both maxima with the primary energy 
is different. Therefore, the consistency between the interpretation of the measurements of X e n ™ x 
and X^ nax in terms of mass composition is not possible for a hadronic interaction model without 
percolation. 



5. Conclusions 

In this paper we have shown that the string percolation phenomenon is expected to appear at 
energies around V* ~ 30 - 40 TeV, i.e, between the LHC energies and those probed in UHECR- 
induced showers collected at the Pierre Auger Observatory. The transition energy depends on the 
assumptions for the proton density profile. 

We have shown that string percolation induces an enhancement in the production of strange 
mesons and the total number of baryons, providing a mechanism to increase the muon content in 
extensive air showers. Moreover, string percolation would harden the particle energy spectrum in 
hadronic interactions and consequently, due to energy conservation, it would reduce the multiplic- 
ity of the highest energy hadronic interactions. 

The effect of string percolation on shower observables was investigated. Percolation leads to 
a reduction of the energy flow into the electromagnetic component of the air shower, thus causing 
an energy shift if the shower should be interpreted with an energy scale with no percolation. The 
natural increment of the energy on the hadronic cascade, caused by percolation, should increase 
the muon content by up to ~ 50% relative to standard proton-induced extensive air showers at 
E = 10 19 eV. 

The depth of the maximum of the electromagnetic shower longitudinal profile displays a tran- 
sition around the percolation energies towards deeper values. The reduction of multiplicity and 
the hardening of the particle spectrum, effectively slows down the growth of the hadronic shower 
causing the depth of muonic component maximum to also become deeper. 
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One of the predictions of the percolation scenario is that the electromagnetic and muon pro- 
duction shower maxima should have different evolutions with the primary energy for a fixed mass 
composition. This could be experimentally checked by combining the information of the Fluores- 
cence Detectors, which measure the electromagnetic with the one provided by the Surface 
Detectors, which can be sensitive to the muon production maximum, X^ nax . 
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